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ABSTRACT 

on" 

Qs. I We show that flows of pure electromagnetic energy are subject to instabilities familiar 

from hydrodynamics and plasma physics — Kelvin-Helmholtz and screw instabilities. In 
i-C ■ the framework of force- free electrodynamics, it is found that electric-like tangential 

["t ■ discontinuities are Kelvin-Helmholtz-unstable. All non-trivial cylindrically symmetrical 

Qv , magnetic configurations are screw-unstable. The Poynting jet of the Goldreich-Julian 

pulsar is screw-unstable if the current density in the polar cap region exceeds the charge 
density times the speed of light. This may be the process that determines the pulsar 
K^ ' luminosity. 

oo 
(N 

^ . Subject headings: instabilities — magnetic fields 

'^ ' 1. Introduction 



o 



It has been suggested that the luminosity of pulsars, Kerr black holes in external magnetic 



c/3 , fields, relativistic accretion disks, and gamma-ray bursts comes originally in the form of a Poynting 

' jet. We show that Poynting jets are subject to Kelvin-Helmholtz and screw instabilities. Generi- 

cally, Poynting jets should lose their stability as they propagate away from the sources because the 
^\ ' guiding magnetic field decreases more quickly than the toroidal magnetic field. Both the guiding 

j^ ■ magnetic field B and the Poynting flux density S decrease as 1/A, where A is the cross section of 

the Poynting jet, and we assume that A grows as the jet propagates away from the source. Since 
the toroidal magnetic fleld and the radial (in local cylindrical coordinates) electric field scale as 
S^'^, these fields decrease slower than the guiding field B. The stabilizing effect of the guiding 
field decreases, causing the jet to lose stability. Instabilities probably lead to a turbulent cascade 
of energy to smaller scales. The ultimate energy sink may be particle acceleration and/or pair 
production. 

We investigate the instabilities in the framework of force- free electrodynamics (FFE). FFE 
has been applied to pulsars, black holes, accretion disks, and gamma-ray bursts (Goldreich &: Ju- 
lian 1969, Scharlemann & Wagoner 1973, Blandford 1976, Blandford & Znajek 1977, Michel 1973, 
Thompson & Blaes 1998, Meszaros & Rees 1997). By definition, FFE applies when the energy- 
momentum tensor is dominated by the electromagnetic fields. Particles carry charges and currents 
but have zero inertia. We present in §2 different formulations of FFE that are used in the stability 



calculations. We show in §3 that electric-like tangential discontinuities are Kelvin-Helmholtz unsta- 
ble. We prove in §4 that all non-trivial cylindrically symmetrical magnetic configurations are screw 
unstable. We give in §5.1 the eigenmode equation for a cylindrical Poynting jet. As an illustration 
of applications of this equation, we show in §5.2 that the Poynting jet of the Goldreich-Julian pulsar 
is screw unstable if the Poynting luminosity exceeds the classical "magneto-dipole" value. 



2. Force-free electrodynamics (FFE) 

FFE is applicable if electromagnetic fields are strong enough to produce pairs and baryon 
contamination is prevented by strong gravitational fields (Blandford & Znajek 1977). Pulsars, 
Kerr black holes in external magnetic fields, relativistic accretion disks, and gamma-ray bursts are 
the astrophysical objects whose luminosity might come originally in a pure electromagnetic form 
describable by FFE. 

FFE is classical electrodynamics supplemented by the force-free condition: 

dtB = -V X E, (1) 

dtE = VxB-3, (2) 

pE + jxB = 0. (3) 

V • B = is the initial condition. The speed of light isc=l;p = V-E and j are the charge and 
current densities multiplied by 47r. The electric field is everywhere perpendicular to the magnetic 
field, E ■ B = 0. The electric field component parallel to the magnetic field should vanish because 
the charges are freely available in FFE. It is also assumed that the electric field is everywhere 
weaker than the magnetic field, E"^ < B^. Then equation (3) means that it is always possible to 
find a local reference frame where the field is a pure magnetic field, and the current is flowing along 
this field. FFE is Lorentz invariant. 

Equation (3) can be written in the form of the Ohm's law. The current perpendicular to the 
local magnetic field can be calculated from equation (3). The parallel current is determined from 
the condition that electric and magnetic fields remain perpendicular during the evolution described 
by the Maxwell equations (1), (2). We thus obtain the following non- linear Ohm's law 

. (B-VxB-E-Vx E)B + (V • E)E x B 

- - — ^ . (4) 



Equations (1), (2), (4) form an evolutionary system (initial condition E • B = is assumed). It 
therefore makes sense to study stability of equilibrium electromagnetic fields in FFE. One can also 
study linear waves and their nonlinear interactions in the framework of FFE (Thompson &; Blaes 
1998). 



If p 7^ 0, E ■ B = follows from (3), if p = 0, this condition is an independent basic equation of FFE. 



It is convenient to introduce a formulation of FFE similar to magnetohydrodynamics (MHD); 
then we can use the familiar techniques of MHD to test stability of magnetic configurations. To 
this end, define a field v = E x B/i?^, which is similar to velocity in MHD. Then E = — v x B and 
equation (1) becomes the "frozen- in" law 

OtB = Vx(vxB). (5) 

From V = E X B/S^, and from equations (l)-(3), one obtains the momentum equation 

at(B2v) = V X B X B + V X E X E+(V-E)E, (6) 

where E = — v x B. Equations (5), (6) are the usual MHD equations except that the density is 
equal B^ and there are order v'^ corrections in the momentum equation. 

3. Tangential discontinuities 

The planar Poynting jet Eq = {V{x),0,0), Bq = {0,U{x),B{x)) is a stationary solution of 
FFE if 

^2 ^ [/2 _ y2 ^ const. (7) 

The eigenmodes of this jet are oc exp{—iujt + ikyy + ik^z) . As shown in Appendix A, the eigenmode 
equation for the x-displacement ^, defined by 6Bx = i(Uky + Bkz)^, is 

{Fa + {u^'-kl-kl)F^ = 0. (8) 

Here the prime denotes the x-derivative, and F = {ujB + kyV)'^ — {k^B + kyll)"^ + {kzV — loU)'^. 

To illustrate applications of the eigenmode equation (8), consider a tangential discontinuity 
Eo = (F,0,0), Bo = {0,U,B) at j; < 0, and Eq = 0, Bq = {0,0, Bq) at x > 0. As follows from 
equation (7), Bq = B"^ + U"^ — V'^. The eigenmode equation (8) is solved by ^ = exp(— k|x|), 
K^ = ky + k'^ — iiP' . From the continuity of F^', F(— 0) + -F(+0) = 0, which gives the dispersion law 

(cjB + kyVf - {k,B + kyUf + {k,V - uUf + (w^ _ ^2)^2 ^ q_ ^g^ 

It follows from the dispersion law (9) that magnetic-like discontinuities, \U\ > \V\, are stable 
and electric-like discontinuities , \V\ > \U\, are unstable Q 



4. Magnetic configurations 

A force-free magnetic field ( V x B x B = 0) with a zero electric field is a stationary FFE solution. 
We will study the stability of the force-free magnetic configurations using the MHD formulation of 



The simplest way to see this is to use a Lorentz boost along z to remove either U or V. 



FFE given by equations (5) and (6). Linear perturbations, denoted b and v, oc exp{—iujt), satisfy 

-iu;b = V X (v xB), (10) 

- iujB\ = V xBxh + V xhxB. (11) 

Define the displacement ^ by v = —iuj$,. From (10), b = V x (^ x B). Now equation (11) can be 
written as 

- uj'^B'^^ = VxBxVx(^xB) + VxVx(^xB)xB = -K^. (12) 

Since the operator K is self-adjoint (Kadomtsev 1966 and references therein), the frequency is 
given by the variational principle 

The equilibrium field B is unstable if the potential energy 

W = f (fr$,K$, (14) 

is negative for some displacements ^. 

Consider the case of cylindrical symmetry, coordinates {r,6,z). The equilibrium field B = 
(0, U{r), B{r)) should satisfy BE' + Ur~^{rU)' = 0, where the prime denotes the r-derivative. For 
an eigenmode oc exp{im9 + ikz), the potential energy reduces to (Kadomtsev 1966) 

W^Jdri/e+ge), (15) 

where ^ is the radial component of the displacement. The other two components of the displacement 

vector, ^g and ^2, were chosen to minimize the energy for a given ^. The functions of radius / and 

g are given by 

{krB+mUf 

g = r-1 AV + m^-l 2 _^piL^^k\^B^ - m^U^)\ . (17) 

^ ^2^2 j^rn? ^ ' (fc2^2 _^ ^2)2 V ) j \ ) 

We now show that magnetic configurations with a non-zero toroidal field U are screw unstable. 
Screw means, e. g., that m = 1, but not m = —1, is unstable for a given sign of k. Kadomtsev 
(1966) gives a clear discussion of the screw instability in plasmas, and shows that the screw mode 
is the most dangerous mode (if the plasma is unstable, it is screw unstable). 

Assume that U{0) = 0, U'{0) > 0, U ^ for r —>■ 00, and U is positive in between. Assume 
that B is everywhere positive. Let k be positive and small. Take m = —1. Let vq be the first 
zero of krB — U. Take a trial function ^ = 1 for r < rQ, and zero at r > tq. We can chose this 
generalized function ^ in such a way that the first term in the energy integral (15) vanishes. The 
second term is an integral of g from to tq. For small k, g = k'^r{krB — U){3krB + U) is negative. 



5. Poynting jets 

Poynting jets become observable only after the instabilities convert the Poynting flux into 
kinetic energy of charged particles. The FFE instabilities may also limit the Poynting luminosity of 
the sources. Thus, the stability of Poynting jets in FFE is of interest for astrophysical applications. 
Unlike the purely magnetic case, the eigenvalue problem is not self-adjoint. We first give the 
eigenmode equation for cylindrical Poynting jets. We then analyze stability of a slowly rotating 
Goldreich- Julian pulsar. We will show that the Poynting jet is screw unstable if the current density 
exceeds the charge density times the speed of light |^. 



5.1. The eigenmode equation 

Cylindrical Poynting jet Eq = (F(r),0,0), Bq = {0,U{r),B{r)), in cylindrical coordinates 
(r, 6, z), is a stationary solution of FFE if 

BB' + Ur-^{rUy -Vr-\rVy = 0, (18) 

the prime denotes the r-derivative. The eigenmodes are oc exp(—iu;t + imO + ikz). Proceeding 
along the same lines as in Appendix A, we obtain the eigenmode equation for the r-displacement 
^, defined by SBj. = i{mU/r + Bk)S,, 

if a = 9i. (19) 



Here 



/ = -rF/K\ (20) 

g = r-^E/K\ (21) 

K^ = fc2 _^ mV~^ - cj^ (22) 

F = c\-cl + ci, (23) 

E={l-K^r'^)F-2 ' """ G-AH-rH', (24) 

G = 2B{loci - kc2) - F, (25) 

H = 4, (26) 

ci = ujB + mV/r, (27) 

C2 = kB + mU/r, (28) 

C3 = kV - ijU. (29) 



^Since charges of either sign are available, the current density can, in principle, exceeds the charge density times 
the speed of light. 



The eigenniode equation (19) reduces to (8) in the planar hmit, r -^ oo, m -^ oo, m/r = 
const = ky. In the purely magnetic case, V = 0, foi to = 0, the functions / and g reduce to (16), 
(17). The eigenniode equation is invariant under the Lorentz boosts in the z-direction - since (to, k), 
{U,V), and (ci,C2) are vectors, one can show that / and g are scalars. 



5.2. Stability of the Goldreich-Julian pulsar model 

To illustrate how equation (19) can be applied, we will derive the stability criterion for a 
slowly rotating Goldreich-Julian (1969) pulsar. The pulsar is called slow if 0,R <C c, where Q is 
the neutron star (NS) rotation frequency and R is the NS radius, c = 1 is the speed of light. The 
pulsar luminosity is carried by a narrow Poynting jet, radius a -^ {Q.R/cY''^R, emanating from the 
polar cap region. Since a <C i?, the jet is approximately cylindrical. In the vicinity of the axis, the 
electrostatic potential is given by y = QBr, B is the NS magnetic field. The toroidal magnetic 
field can be written in the form U = 0,u{r)Br. We will assume that ^^[/(r) decreases from some 
positive value 0[/ at r = to at r > a. The value of Of/ is not known, but it was suggested that 
Qc/~0. 

The following analytical derivation of the sufficient condition for screw instability has been 
confirmed by numerical simulations of the eigenniode equation. Take A; = (1 — e)Qu, where e is a 
small positive number. If fi = 0, the screw mode {m = — 1) is unstable and has a small growth rate. 
What happens when we "turn on" Q1 Since f/g ~ 1/fc^ ^ a^, the eigenniode is approximately a 
step- function, ^ = 1 at r < vq, ^ = at r > rQ, where tq is the first zero of C2, and the real part 
of the eigenfrequency is Re w = —17, so that |ci| <^ |c2|. The instability criterion is similar to the 
purely magnetic case of §4, 

/ drg < 0. (30) 

Jo 

In the leading order in QR/c, 

g={k^ - n^){kB - U/rf + 2{k^ - n^){k^B^ - U^/r^) - 3n^{kB - U/rf. (31) 

The second term in g is the largest in absolute value. This term is negative if fi < Q.u. Therefore 
the jet is unstable so long as fi < Viij. It follows that the luminosity of a stable Goldreich-Julian 
pulsar does not exceed L ~ a^cUV ^ {U,R/c)'^cB'^R? - the "magneto-dipole" value. 

I thank John Bahcall and Roger Blandford for discussions. This work was supported by NSF 
PHY-9513835. 



A. Planar symmetry 

We derive in this Appendix the eigenmode equation (8). Consider a planar Poynting jet 
Eo = {V{x), 0, 0), Bo = (0, U{x), B{x)), BB' + UU' - VV = 0, the prime denotes the x-derivative. 



The eigenmodes are oc exp{—iujt + ikyy + ikzz) . The FFE equations (l)-(3) for the perturbed fields 
and currents (denoted e, b, j) are 

UJOx = KyGz ~ kzGy, \-^^) 

ujby = kzCx + ie'z, (A2) 

Ldbz = -kyCa: - le'y, (A3) 

Loe^ = -kybz + kzby - ij^, (A4) 

ujcy = -kzbx - ib'z - ijy, (A5) 

ojez = kybx + iby-ijz, (A6) 

Vex + V{e'x + ikyCy + ifc.e,) - B'bz - U'by + Bjy - Ujz = 0, (A7) 

V'ey + U'bx - Bjx = 0, (A8) 

V'cz + B'bx + Ujx = 0. (A9) 



From (A8), (A9), 
From (A4), (A8), 
From (A5), (A6), (A7) 



Vbx + Uey + Bez = 0. (AlO) 

V'cy + U'bx = iB{u!ex + kybz - kzby). (Alf) 

{Bbz + Uby - VcxY = i{ciey + cs^^ + cge^), (A12) 



where we denoted ci = toB + kyV , C2 = kzB + fc^f^, C3 = kzV — toU . For further derivations the 
following relationships are useful: ciU — C2V + C3B = 0, and cikz — C20J — c^ky = 0. 

We solve (Al), (AlO) by introducing the x-displacement ^ (called the x-displacement because 
according to (A13) bx is obtained from U and B by an x-displacement S,)'- 

bx = ic2Ci dy = -icii, dz = -ic3^- (A13) 

Define 

r, = Bb^ + Uby - Vex. (A14) 

Equation (A12) takes the form 

v' = {cl-4 + 4)^. (A15) 

From (A2), (A3), (A13), (A14), (All) 

B{ciS,y = -{loV + kyB)ex + ujUby - tor], (A16) 

{csO' = -kzex+ujby, (A17) 

{C2U' - ciV')^ = c\ex - C2by + ky-q. (A18) 



Excluding e^ and hy from (A16)-(A18), we obtain 

(c? - ci + ci)^' = ikl + A:2 - u?)^. (A19) 



From (A15), (A19), 



((C? - C^ + Cl)^')' = {kl + fc2 _ ^2)(^2 _ ^2 ^ ^2)^^ (A20) 
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